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^-H Abstract. - Non-Markovian processes have recently become a central topic in the study of open quantum 

systems. We realize experimentally non-Markovian decoherence processes of single photons by combining 
time delay and evolution in a polarization-maintaining optical fiber. The experiment allows the identification 
' ^ ' of the process with strongest memory eff'ects as well as the determination of a recently proposed measure for 

^ the degree of quantum non-Markovianity based on the exchange of information between the open system and 

I its environment. Our results show that an experimental quantification of memory in quantum processes is 

^ indeed feasible which could be useful in the development of quantum memory and communication devices. 



Introduction. - The theory of open quantum systems de- 
scribes how a system of interest is influenced by the interaction 

^sf-jwith its environment ||T|. This interaction often leads to a loss 
of the quantum features of physical states and has a great im- 

CN pact on the dynamical behavior of the open system due to the 
non-unitary character of the time evolution. Since any realistic 
physical system is coupled to its surroundings, open quantum 
systems and their description plays an important role in many 
applications of modern quantum physics. 
^ While Markovian, or memoryless, quantum processes are 

k^well understood in the framework of the theory of quantum 
Markovian master equations developed during the 70's |[2j[3|, 
non-Markovian processes with memory have recently become 
of central importance in the study of open systems |[4}|T2|. In 
general, non-Markovian features can arise, e.g., because of a 
strong system-environment interaction, structured reservoirs, 
due to initial system-environment correlations, or couplings to 
a low- temperature environment or spin bath. Moreover, recent 
developments in experimental technology allow reservoir engi- 
neering 1 13], study of quantum correlations tl4J , and the de- 
velopment of quantum simulators for open systems (TSj . 
On one hand, non-Markovian systems are not yet well- 
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understood even on the fundamental level, and many eflTorts 
have been devoted to the development of theoretical tools for 
the treatment of quantum dynamics with memory |[6}fT2|. As a 
matter of fact, the very definition of quantum non-Markovianity 
has been recently under a vivid discussion |T0}|T2| . On the 
other hand, there are indications that non-Markovianity may 
play a role, e.g., in energy transport in certain photo synthetic 
complexes |[T6[[T7 1 and can be exploited for quantum metrol- 
ogy fT8| and cryptography 1 19 |. 

In this Letter we implement non-Markovian processes for 
single photons and demonstrate that it is possible to determine 
experimentally the amount of memory in the system. Our work 
is based on a recent theoretical proposal to study the informa- 
tion flow between the system and its environment, and in par- 
ticular to quantify non-Markovianity with the help of the back- 
flow of information from the environment to the system | [TT| . 
This allows us to identify a specific process having the largest 
memory among all experimentally implemented processes, and 
then to quantify the amount of non-Markovianity for this pro- 
cess. There exist earlier experimental works on dephasing in 
photonic systems pO)|2T|, and results which may be regarded 
as indications for non-Markovian behavior |22,23|. However, 
the recent theoretical progress is opening the path for rigorous 
experimental tests of quantum non-Markovianity, for example, 
by controlling the initial state of the environment p4| , or by 
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modifying the interaction between the system and the environ- 
ment, as is done here. We also note that a theoretical proposal 
to witness initial system-environment correlations by studying 
information flow between the system and the environment |25j 
has been recently realized experimentally |26, 27|. 

Theoretical framework. - We consider a pure dephasing 
quantum process for which the density matrix p of the open 
system evolves according to the master equation 

^ = -/f K,P«].f (c..P«c..-p(0). (1) 

Here, €(t) represents the time-dependent energy shift and y(t) 
the time-dependent rate of the decay channel described by the 
Pauli operator cr^. Generally, non-Markovian dynamics can be 
described, e.g., by memory-kernel equations or time local mas- 
ter equations where the decay rates depend on time ||1J, and we 
have chosen the latter due to its conceptual simplicity for the 
current purpose. Recently, time local equations have received 
a great deal of attention |[6][8][TT][T2][28j, and for a single chan- 
nel system, as is the case here, one can directly associate the 
non-Markovianity with the appearance of negative periods of 
the decay rate |8|pl|. 

The open two-state system we consider consists of the hori- 
zontal and the vertical polarization states of a photon, \H} and 
|y), respectively. The dephasing process influences the co- 
herences between the polarization components and the evo- 
lution can be described by the decoherence function K(t) = 
exp[-^ dfiy{f) + ie{f))\ which is connected to the energy 
shift and the decay rate of the master equation ([T]) by the rela- 
tions 
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y{t) = -%[k(t)im. (2) 



The corresponding dynamical map which maps the initial 
polarization state p(0) to the state p{t) = O^p(O) at time t is then 
given by 

Ph,h(0 = Ph,h(0), Pv,v(0 = Pv,v(0), 
Ph,v(0 = K*(t)pny(0), Pv,h(0 = /^(OPv,h(0). 

The environment consists of the frequency degrees of free- 
dom \<jij) of the photon and we consider an initial product state 
between the system and the environment Ptot(O) = P(0) ^ 
JdojJ doj'A(oj)A''(aj')\aj}{oj'\. The function G(aj) = \A(aj)\^ 
represents the normalized frequency distribution of the pho- 
ton which in our experiments is a Lorentzian distribution with 
central frequency ojq and full width at half maximum 26co 

(FWHM), G(oj) = ^- \r^' We investigate the case in 

which the evolution of the frequency distribution depends on 
the polarization state such that the total system dynamics is de- 
scribed by the unitary operator 

+ |V)<y|0 J doje-'''''''^'^'\oj}{oj\, (3) 
where the factors uu(t) and wv(0 may depend on time. 
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Fig. 1: (Color online) The experimental setup consisting of three 
parts: Preparation of a single photon state, evolution in the delay setup 
followed by the PM optical fiber, and state tomography (Tl and T2). 
Used abbreviations: BS, beamsplitter; LP, long pass filter; IF, inter- 
ference filter; HWP, half- wave plate; PBS, polarization beam split- 
ter; QWP, quarter-wave plate; PM, polarization maintaining; APD, 
avalanche photodiode. 



In the experiment we determine the measure for the degree of 
non-Markovianity constructed in Ref. |TT| . For a given quan- 
tum process of an open system described by a dynamical map 

this measure is defined by 



7V(0) = max ( 



dt cr(r,pi,2(0)). 



(4) 
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Here, p 1,2(0) are two initial states of the open system and 
<^(^P 1,2(0)) = j^D(pi(t),p2(t)) is the rate of change of the 
trace distance D(pi,p2) = |tr|pi - P2I. The trace distance 
represents a measure for the distinguishability of two quantum 
states. In Eq. ^ the time integral is extended over all inter- 
vals in which the trace distance increases and the maximum is 
taken over all pairs of initial states. The measure A/'(0) thus 
quantifies the maximal total increase of the trace distance dur- 
ing the time evolution, which can be interpreted as the maxi- 
mal total amount of information that flows from the environ- 
ment back to the open system fTT| . In our experiment the in- 
crease of the trace distance signifying non-Markovian behav- 
ior is restricted to a single time interval [toJi] which yields 
Nm = maxp^2(0) [^(Pi(^i),P2(^i)) - D(pi(toXp2(to))]. 

Experimental framework. - We measure the non- 
Markovianity of a quantum evolution process given by the mas- 
ter equation ([T]) which is implemented for single photons emit- 
ted from a quantum dot (QD). The experimental setup is shown 
in Fig. [T] and includes three parts: Preparation, time evolution, 
and state tomography. The first part consists of the generation 
of a single photon and its preparation in the pure initial state 



1^) = cos(0) \H} -h sin(0) \V} , 



(5) 



where is the angle of the polarization direction of the pho- 
ton from the horizontal direction. A self-assembled InAs/GaAs 
QD sample (at a temperature of 7K) [29 , 30] provides the sin- 
gle photon source. A He-Ne laser pumps the sample through a 
beam splitter (BS, high transmission efiftciency of 92%) and a 
50X objective. The photon is collected by the same objective 
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Fig. 2: (Color online) (a) Difference of the trace distances at times tf 
and to for various delay lengths xq. The initial states are prepared with 
a fixed pair of angles (135^, 45^). (b) Final trace distance at time tf for 
the pair (135^, 45^). The red solid line is the theoretical curve (see the 
text) giving l/Soj = 35.8 + 1.9ps. 

and separated by a 785nm long pass filter (LP) and an inter- 
ference filter (IF, the central wavelength of the pass band is 
tunable around 950nm and the FWHM is 0.7nm). Then, a po- 
larizer and a half- wave plate (HWPl) prepares the photon in 
the initial state ([5]), where can be changed by rotating HWPl. 

The time evolution consists of two main contributions: A 
delay setup and a polarization maintaining (PM) optical fiber 
with a HWP3 between the two contributions. The delay setup 
includes a polarization beam splitter (PBSl), two 45^ (from 
the horizontal direction) placed quarter- wave plates (QWPl 
and QWP2), and two mirrors whose relative delay can be con- 
trolled. The decoherence function can generally be written as 
K(t) = J dojG(oj)Qxp(iojt). Considering the delay setup and 
using the earlier mentioned Lorentzian photon frequency dis- 
tribution, we find 



K(t) = eXp[/6L)o^ - 



(6) 



where we use the notation t = 2x/c with x corresponding to the 
delay length between the mirrors. The maximal value of the 
delay for a given set of experimental runs is denoted by xq and 
the corresponding time by to = 2xo/c. 

In the second phase of the evolution for times t ^ t^ the 
photon travels in an optical fiber with length / = 100m and a 
designed birefringence of A^o = 3.5 x 10"^. The total sys- 
tem evolves according to the Eq. ^ where now unit) = riu/n 
and uyit) = fly In, with h denoting the average of the ordinary 
and the extraordinary refractive indices nw and ny, respectively. 
The birefringence of the PM fiber influences the relative phases 
of the polarization components I^O^ 21J and by tracing out the 
environment leads to the decoherence function 



Kit) - exp 



icoo to 



An 



(^- ^o) -^^^ko 



An 



-(t-to)\ 



(7) 



which is used for times ^ ^ ^ , where tf is the termination 
time of the process when the photon exits the fiber, and An is 
the real birefringence. 

In the third and last part of the setup a tomography [ 31 1 of 
the final states, i.e., of the states at the end of the fiber is car- 
ried out (Tl). We also perform a tomography of the states after 
the delay setup (T2). Both tomography setups are identical uti- 
lizing two single-photon avalanche photodiodes (APDs). The 
count rate is about 7000/s, and the integration time is 4s. 
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Fig. 3: (Color online) Difference of the final and initial trace dis- 
tances. The initial states are prepared with a fixed delay xq = 19. 15mm 
and various pairs of angles iO,^). A maximum value is observed for 
i6, ^) = (135^, 45^). The inset shows an elaborate measurement around 
the pair(135^45^). 



While the length of the fiber is fixed, we can control the max- 
imal delay to before the photons enter the fiber. The evolution 
of the reduced system follows the master equation ([T]) for all 
times, where yit) is given by 



r(f) 



Soj, 

n ' 



^t <to, 
to ^ t < tl, 
tl ^t < tf. 



(8) 



These can be obtained by using Eqs. ^ and ([6][7]), and ti = 
(1 -h n/An)to > to. We can control the length of the negative 
decay period by changing ^o- Note that if to is small enough, 
then tf < tl and the last positive period does not appear since the 
decay rate is negative during all of the fiber evolution. It is also 
easy to see from Eq. (|7| that the period of a negative decay rate 
corresponds to an increase of |a:(0|. By fixing to, corresponding 
to a fixed delay xq, one also fixes the dynamical map which 
we denote by ^ ^nd there exists a diflferent process for each 
value of xo. We are now interested which one of the processes 
O^o has the highest value of non-Markovianity and what is the 
corresponding value of xq. 

Results. - To identify the process having the largest value 
of non-Markovianity, we first fix the angles of the initial- state 
pair as (6>,^) = (135^45^) [see Eqs. Q and Q]. After the 
delay time to, the photons travel in the fiber and a state tomog- 
raphy is performed both before and after the evolution in the 
fiber. The diflTerence of the trace distances at the correspond- 
ing times tf and is denoted by ADiO, ^, xo) and the result is 
shown in Fig. [2ja). We find that the maximum is located at 
Xo = 19.15mm and takes the value AZ)(135^45^ 19.15mm) = 
0.962 ±0.01 1 . If we subtract the dark and background counts of 
the APDs (about 150/s), this value is corrected to 0.988+0.011. 
The uncertainty is due to the counting statistics. Fig.[2jb) shows 
the final trace distance Dt^ as the function of the used delay 
Xo. By fitting the experimental data to the theoretical curve 
Af (-^o) = exp(-^a;|A^/ - 2xo\/c) [the red solid line in Fig.[2jb)] 
we determine the parameter l/Soj = 35.8 ± 1.9ps. 

For Xo ^ 19.15mm the decay rate is positive when t < 
to and negative for all times the photon travels in the fiber, 
^0 ^ ^ ^ = ^1 . Hence, the results in Fig. |2] (a) show that 
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Fig. 4: (Color online) (a) The measured spectrum of the quantum 
dot. Fit with Lorentzian curve gives the width 1 /Soj = 34.0 + 1.7ps.(b) 
Simulation results for the quantity AD(6,^, 19.15mm) corresponding 
to the experimental data of Fig. [3] The colors represent the values of 
AD. 



for xo ^ 19.15mm, the process o^^-^^"^"^ gives the maximal 
non-Markovianity within this group of processes o^^^^^-^^"^"^. 
On the other hand, for the processes with xq ^ 19.15mm, 
the smaller xq the shorter is the negative decoherence period 
since the negative decay occurs for ^ ^ ^< < ^f- Con- 
sequently, the point xq = 19.15mm also gives the maximal 
non-Markovianity for values xq ^ 19.15mm. Thus, we can 
conclude that (^^^-^^^^ yields the maximal non-Markovianity 
among all processes experimentally implemented. 

In the second step, we fix the delay corresponding to 
the process O^^-^^"^"^ and carry out the maximization over 
pairs of initial states in the definition ^ of the measure by 
changing the pair of angles (0,^). Experimental results for 
AD(6,^, 19.15mm) are shown in Fig. [3] We clearly see that 
AZ)(135^, 45*^, 19.15mm) still represents the maximal value 
among all pairs of initial states of the form given by Eq. ([5]). We 
do not perform the maximization over the whole state space, 
because restricting to pairs given by Eq. ([5]) is suflftcient for 
finding the maximizing pair as we will see later on. 

Let us discuss this point in more detail. Firstly, we find 
that 2 X 19.15mm is approximately equal to A^o^ = 35mm. 
The discrepancy is explained by noting that the birefringence 
given by the manufacturer of the fiber at its designed wave- 
length 780nm is A^o = 3.5 x 10""^. However, the wavelength of 
our chosen peak is at 946. 3nm and consequently the real bire- 
fringence in our experiment is A^ = 3.83 x 10"^ [the value 
used in Eq. ([7])]. This suggests that the birefringence eff'ect of 
the PM optical fiber compensates exactly the dephasing caused 
by the delay setup when xq = 19.15mm. This means that the 
open system has recovered all the information that it lost dur- 
ing the earlier part of the evolution. The experimentally deter- 
mined value for the measure of non-Markovinity is given by 
MO) = 0.962 ± 0.012, and by 7V(0) = 0.988 ± 0.01 1 when ac- 
counting for dark and background counts. This result matches 
very closely the theoretical result 7V(0) = 0.972 obtained by 
using Eqs. ^Sff\ and the above initial states. The diff'erence be- 
tween the experimental and theoretical result may be caused by 
mode flips and temperature changes in the PM fiber. Notice, 
that increasing the fiber length and correspondingly the delay 



time ultimately give an increase of 1 for the pair of angles 
(135^,45^). Since the trace distance is by definition bounded 
between and 1 there exists no pair exceeding this value of 
increase and thus the pair with the angles (135^,45^) must be 
maximizing. 

We can also directly measure the spectrum of the quantum 
dot, see Fig. |4ja). The obtained width l/6co = 34.0 ± 1.7ps 
matches very well the theoretical value I /Soj = 35.8 ± 1.9ps 
[see Fig.[2jb)]. Moreover, we have performed numerical sim- 
ulations of AD(6,^, 19.15mm) shown in Fig. |4jb). The found 
maximum value at (135°, 45°) coincides with our experimental 
result in Fig. [3] 

Conclusions. - We have experimentally realized a fam- 
ily of pure decoherence processes in a non-Markovian open 
quantum system. The setup allows to identify the process with 
strongest memory eff'ects in this family and to determine the 
corresponding measure for non-Markovianity (|4]). This clearly 
demonstrates the measurability of this quantity, including in 
particular the maximization over the initial states. Our results 
show that wide scale experimental studies on fundamental and 
practical aspects of quantum systems with memory are indeed 
becoming feasible in the wake of the recent vivid theoretical 
debate. Moreover, our results may be helpful in the study of 
quantum communication processes [32J . 
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